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height of spine (m) 

specific heat of surrounding fluid (J/kg: Kk) 
diameter of spine (n) 

differential amount cf heat (w) 

cross sectional area function for spine (pm ) 
profile function for spine (nm) 

perimeter function fcr spine ({m) 
gravitational acceleration (n/s*) 
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heat transfer coefficient (w/m + K) 

thermal conductivity (w/m-K) 

thermal cond. of surrounding fluid (w/m-K) 
thermal cond. of insulation (w/m- K) 
thermal cond. of spine material (w/m-K) 
Tin heat transfer parameter (1/n) 

Nusselt number (dimensionless) 

Prandtl number (dimensionless) 

heat flow into spine at base (w) 

Rayleigh number (dimensionless) 

Reynolds number (dimensionless) 
temperature of spine (fin) surface (K) 
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top Tw 


Z (2) 


temperature of surrounding fluid (K) 


enhancement factor a function of the roll angle 


GREEK SYuBers 


a combination of variables defined by egn (A.9) 
volumetric coefficient of thermal expansion (1/K) 
diameter of spine (m) 

roll angle (degrees) 

temperature excess ( k) 

temperature excess at spine base ( K) 

temperature excess at spine tip ( K) 

dynamic viscosity of surrounding fluid (kg/m-s) 


kinematic viscosity cf surrounding fluid (n*/s) 
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I. INTRODUCTION 

Electronicaily steered phased array antennae frequently 
exhibit degraded performance due to the lack of proper heat 
transfer between the interior crystalline materials and the 
environment. This causes higher than optimum temperature 
environments for the crystalline materials and a subsequent 
loss of reliability[ 1]. 

One possible solution to such a problem considers a 
surface augmentation through the use of pin fins (or Spines) 
to assist in the dissipation of the heat. 

The problem 1S more severe when the installation is 
Stationary so that the dissipation from the spines is in the 
natural convection mode and, because of this, the focus of 
this investigation is in the area of natural convection heat 
transfer and methods of enhancenent in this domain. 

Convective heat transfer from horizontal cylinders has 
been the subject of extensive studies which were summarized 
Dy sion Gans. as In addition, studies of heat exchanger tube 
bundles are summarized by Incropera and and DewWitt[3]. 
However, after a considerable search of the literature, only 
a small amount of experimental data and analytical correla- 
tions dealing with natural convection in an array of hori- 
zontal cylinders, has been found. 

Intuitively, one feels that the orientation of the array 
must have an effect on the heat transfer but apparently no 
known experimental correlations have been formulated. fThus, 
the focus of this thesis 1S in the effect of roll angle on 
the heat transfer rates from an array of cylindrical spines. 

Additional impetus for investigations in Similar areas 
with cylindrical Spine arrays derives from the electronics 
industry where such cooling arrangementS are used in the 


external cooling of dense integrated circuit arrangements. 


10 


Dimensionless correlations using dimensionless parame- 
ters are usually developed for presentation of experimental 
results and this allows for pcssible application to other 
systems. In develofing the solution for the heat transfer 
across a boundary layer, nets has been shown by 
Bentachting(/4}, and others [3,5], that the Nusselt number 
(Nu) iS a proper dimensionless parameter for natural or 
forced convection. In forced convection, the Nusselt number 
is a function of the Reynolds number (Re) and the Prandtl 
number (Pr). However, in pure natural convection where the 
buoyancy is the only driving fcrce, the fluid velocity is 
determined entirely by the quantities contained in the 
Grashof number (Gr) which is used instead of the Reynolds 
number to establish dynamic similarity. The Grashof number 
is developed guite thoroughly by Kreith and others [3,5] and 
guickly appears in a dimensional analysis. 

In addition, for gases the Frandtl number (Pr) is nearly 
constant so that the Nusselt Modulus (Nujy) can be considered 
to be a function of the product of the Grashof number and 
the Prandtl number. This product is defined as the Rayleigh 
or number (Ra) and is discussed to further detail in the 
following paragraphs. 

In order to establish the heat transfer coefficient and 
to develop a correlation of the data with dimensionless 
parameters, relationships for the heat transfer coefficient 
and the Nusselt number were eStablished and measurements 
were taken which permitted its calculation. A summary of 
the derivation is provided in Figure 1.1 anda detailed 
derivation is provided in Appendix A. The relation holds 
for the condition that no heat transfer occurs at the fin 
ED. Such a condition can be effectively achieved by insu- 
Hating the tip. 

After calculating the heat transfer coefficient the 


Nusselt number was then determined from the relation: 
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Pigure 1.1 A Relation for the Heat Transfer Coefficient 
‘ ave 
ey a Qe 


where Nu, is the Nusselt number based on the spine dlamemeans 
The Rayleigh nunbec was deternined from the product of 
the Grashof nuuwber and the Prandtl number which were calcu- 


lated from the tollowine: 


5 BGuery acs 





Gr = 
v2 


Ra = GrPr Clie-2) 


where again the subscript "d" infers a dependence on the 
Spine diameter and where the symbols are as defined in the 


nomenclature list. 
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If. EXPERIMENTAL APPARATUS 
Figure 2.1 displays the configuration of the experi- 
Mental apparatus and the technigues used in data acguisi- 
Ean. The following paragraphs provide a more indepth 


description of the component parts. 


A. TEST DISK 


The disk, shown in Figures 2.2 and 2.3, was fabricated 
entirely of 2024-T4Y aluminum plate,6.35 cm in thickness. A 
Circular disk of diameter equal to 22.86 cm was cut from 
this plate and fitted with an array of 16 cylindrical sfines 
as the heat transfer surfaces. Six equally spaced radial 
holes were drilled around the perimeter of the disk for 
insertion of six 200 watt cartridge heaters which provided 
the heat input necessary to elevate the temperature of the 
Spines. The heaters were powered by a controllable power 
source in order to attain a steady-state condition for data 
measurements. The disk was fully insulated with soft rubber 
refrigeration insulation (k,; = 0.15 w/m-K) to minimize heat 
losses from the disk surface due to both radiation and 


CONVECEIOR. 


BH. CYLINDRICAL SPER ES 


The spines, Shown in Figures 2.3 and 2.4 were similarly 
constructed of the 2024-T4 aluftinun. Two sizes of spines 
were used in data gathering (0.9525 cm diameter x 3.75 cm 
length and 1.27 cm diameter x 5.08 cm length) to provide a 
wider range of Rayleigh Numbers. Extension pieces, also 
made of the same aluminum, were necessary to extend the 


actual base of the heat transfer surface (the cylindrical 
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Figure 2.2 Test Disk and Spine Assembly 


spine) to the outer level (exterior surface) of the 


InStiLa tom. 


C. THERBOCOUPLES 


All thermocouples were fabricated from copper-constantin 
(Type-T) of AWG No. 30 wire. They were butt welded witha 
spherical bead of approximately 0.05 cm diameter and 
staller. All interconnecting cables were also of copper- 


constantin in order to minimize the uncertainty and to 


preclude, as much aS possible, spurious thermoelectric 
Cieects:. All thermocouples were connected to the ACUREX 
Autodata Nine data acyuisiticn system for temperature 


EGCOL danige 
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Figure 2.3 Sketch of Front View of Test Disk 


D. ACUREX AUTODATA NINE 


The thermocouple modules in the data acguisition 
system's thermocouple modules include an ice-point compensa- 
tion network and _= frovide automatic data scanning and 
Peentout for up to 100 channels of input. Due to the 
internal ice-point ccmpensation networks, no external ice 
baths or reference junctions were required. Pn place 7cala— 
bration checks were made at various temperatures over the 
expected temperature range in crder to develop temperature 
Peeeection factors for each thermocouple, its connecting 
Seeing and the Autodata Nine circuitry as a unit. in addi- 
tion, the complete disk and spine array was immersed in an 


ice-Lbath in order to establish cne of the check points. 
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E. ENVIRONMENTAL CHAMBER 


The test system was housed in an environmental chamber 
in which the ambient temperatures could be maintained over a 


range between 35 and 135 degrees Fahrenheit. 


A. DATA GATHERING 


The experimental procedure was to set the disk at the 
desired roll angle using a liguid bubble level indicator and 
protractcr to measure the angle. The temperature control 
thermostat of the environmental chamber was then set to the 
desired temperature. Steady state environmental conditions 
were attained in approximately six to eight hours. Next the 
heaters imbedded in the disk were energized to elevate the 
temperature of the disk and base of the spine to a tempera- 
ture chosen for a predetermined Rayleigh number. The input 
voltage for the heaters waS varied and the temperatures 
monitored until a steady state condition was reached. Ths 
waS indicated by constant input power and constant tempera- 
tures. The time to attain a total steady-state condition 
waS approximately three to five hours in addition to the 
time reguired for the chamber tc stabilize. 

A steady state condition was assumed to exist when no 
change in temperature occurred for four consecutive tempera- 
ture readings taken at fifteen minute intervals. These 
temperatures were then used for the input data for the 
calculations. Typical dataruns are presented in Appendix B. 

After completion of data taking at a set condition, 
changes were made in the roll angle, the temperature of the 
system or temperature of the environment, and the process of 
attaining and maintaining a new condition of steady state 


waS again started. 
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B. DATA REDUCTION 


Temperature data printed by the Autodata Nine during the 
individual runs was used aS input parameters (corrected fron 
the calibration curve for each thermocouple) for calculation 
of the desired parameters for plotting and evaluation. 

All desired output parameters were calculated by the 
computer using appropriate relationships for the thermophy- 
Sical properties [3] as functicns of temperature and meas- 
ured temperatures aS input paraueters. The computer program 
used for the calculations is provided in Appendix E. 

Graphical display of the data was also utilized to aid 
in data analyses. These curves are shown aS Figures B.1 
ERE OUGH Bs -« 
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IV. RESULTS AND CORRELATION 


= Se = ao ——— ee ee ee 


Ae RESULTS 


During each data collection run the temperatures at the 
Base and tip of cach spine, as well as four environmental 
temperatures were recorded. These temperatures were used as 
inputs for calculation of the various output parameters for 
each spine. These sixteen calculations were then averaged 
to provide the data feoint for that set of conditions. Forty 
seven data points were obtained. Typical data runs are 
included as Appendix C. 

At a specified angle, the natural logarithm of the 
Nusselt number was flotted as a function of the natural 
logarithm of the Rayleigh number for each of the data pcints 
at a specified angle. These flots are included as Figures 
B.1 through B.5 in Appendix B. This graphical display was 
then used in order to facilitate correlation of the test 
ote « 

The plotted results indicated that a small dependence on 
roll angle did exist. Figures B.6 and B.7 were plotted to 
Show these effects more specifically. Figure B.7 was 
plotted only for Rayleigh numbers between 579 and 586 
because of the narrowness of the overall Rayleigh number 


range obtained. 


B. CORRELATION 


For each of the curves showing log (Nu) as a function of 
log (Ra), a least squares curve fit routine was used anda 
linear relation established tc approximate these curves 
following a method presented by Kraus and Bar-Cohen[1]. A 
correlating equation was obtained for each of the froll 


angles (%), where Z(@) defines a function of the roll angle. 
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The values ot 2(f) are tabulated in Table le 


TABLE 
Enhancement Factor 
evo) EAL Enhancement 
Angle(@) Factor (€Z60)> 
Oh Bee5 9 6 
POS 0.5986 
202 O607 5 
SO 0254 50 
BS ee 0.5928 


The values of Nusselt numbers obtained from the experi- 
mental results are a functiot of a given Ravleiyh nuwber and 


roll angle and way Ke correlated via: 


O22 36 
Nu = Z(@)Ra (4.1) 


However, orserve that the dependence on the roll angle is 
weak (Table 1) and no eathenmatical relationship was deter- 
mined for the values or Z(f). 

Results were also compared with several previously 
published correlations for the single horizontal cylinder 
Churchill and Chu{6] suggested an ejuation Of tre toma 


Z 


Orel eee 


(9/16) 7] (8/27) 
, (Pres 





ebe, 


ae 


which correlated experimental results of Several researchers 
very well. | 

The values obtained in this experimental program agreed 
favorabiy with these predicticns with an average of 19% 
deviation throughtout the range investigated. Noting that 
the effect of heat transfer by radiation were also present, 
the deviation is attributed to radiation losses primarily. 
These losses would aprear to nct have a detrimental effect 


on the observed roll angle dependence. 


Zs 


A. CONCLUSION 


The results of the experimental program indicated that a 
dependence on roli angle did exist. This dependence is 
relatively minor (i.e., on the order of 3%), however for 
temperatures in the range of z00 degrees Fahrenheit this 
could mean a dirfference of 5 to 6 degrees in the surface 
temperature of an operating system. Therefore, further 
investigation at higher Raylé€igh numbers nay be deemed 


appropriate. 


Be. RECOMMENDATIONS 


The following reccmmendations are provided for possible 


follow on projects of a Similar nature. 


1. Conduct the experimental prcgram with the goal of deter- 


mining and removing the radiation heat transfer. 


2. investigate other devices and methods of temperature 
measurement in order to remove more of the uncertainty 


associated with the use of thermocouples. 
3. Conduct a Similar experiment for the pitch angle. 


4. Attempt a Similar program with water as the convecting 


medium. 
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APPENDIX A 
DEFINITION OF THE HEAT TRANSFER COEFFICIENT 


Fins and spines (pin fins) of various geometries and 
thermal conductivities respond differently to identical and 
uniform heat sources and _ sinks. Similarly, there are 
numerous wayS in which the temperatures and heat transfer 
coefficients of sources and SixrksS may vary. Important to 
the analysis of fin geometries are the constraints or 
assumptions which are employed to define and limit’ the 
probiem and often to simplify its solution. The analysis of 
fins and spines emfloy the assumptions of Murray{7] and 
Gardner[8]. These limiting assumptions are: 

1. The heat flow in the fin and its temperatures remain 


constant with time. 


2. The fin material is homogeneous, its thermal conduc- 
tivity is the same in all directions and it remains 


constant. 


3. The heat transfer coefficient between the fin and 
the surrounding medium is constant and uniform over 


the entire surface of the fin (spine). 


4. The temperature of the medium surrounding the fin is 


uniforn. 


5. The fin thickness 1S so shall compared with its height 
that temperature gradients across the fin thickness 


may be neglected. 
6. The temperature at the base of the fin is uniforn. 


- There iS no contact resistance where the base of the 


fin joins the prime surface. 


25 


Shc There are no heat sources within the fin itself. 


9. The heat transferred thrcugh the outermost edge of 
the fin 1S negligible compared with that leaving the 


fin through its Lateral supiacc. 


10. Heat transfer to or from the fin iS proOpOortiond lees 
the temperature excess between the fin and the sur- 


rounding mediup. 


In addition to proposing a frofile function for longitu- 
dinal fins Gardner[{8] proposed a profile fuUnct1On aes 


spines: 





f(x) = = (#F" ‘ (A-1) 


With the appropriate value of n, eguation A.1t may be used 
for the development of a generalized differential eyuation 
for spines. 

Figure A.1 showS a Spine of arbitrary profile. Lt Cam 
be seen that the spine cross section normal to the flow of 
heat, the confining profile, and the perimeter of the spine 
are perfectly arbitrary functiors of the distance x £rommuae 
tip of the spine. A differential equation for the tempera- 
ture excess, O=t-t, may te written by considering the 
heat conduction into and out of tke element dx through the 


cross section £ (x): 


Zc 


<< 
it 
o- 
oc 

il 

© 


i ee ce ec et ee sc ee ce em eee me ee ee eee 


Pigure A.1 Spine of Arbitrary Profile 


The difference in heat flow into and out of the element 
must equal that dissipated by the surface of the spine. iE 
Meostipation occurs ty convection and h is the convection 


@oetricient 


igi hh f(x) Q dx 


where £4 (x) defines a perimeter function P(x) which depends 
Upon the distance x fron tho oriyin of the coordinate 


Systen. When convection and conduction are e‘juated, 


d dQ} _ 


Upon reatrangenent the generalized ‘itferential equation 


becomes 





ee az9 df, (x) do - 
x) ye + a) 
al x 
dx dx dx k 3 (Ae 2) 


The celatlonship heeweon es tame f, (4) LS 


fox) = 47, [F(x] 


Where £,(x) was defined by equation A. 1. Thais @Guat 1 Onegai 


Muy ce written 


[f,¢x)]° a*o d f.(x}*de 2h FA (x)8 | (A. 3) 
Scag ta es —_ — —_ = 0 
dx dx dx k 


Wiere 


Equation A.3 1S a4 Second-order ditferential eG wae 
With variaule cceffiicieLtS except when the Spine “@mgee 
Section notmal to Neat flow is “comctani. it aay he solved 
by termwise comparison with the generalized Bessel equataioge 


For tke cylinirical Stine, shown in Figure 
d 
Pe. 
Z 2 


and 





d 2 


Por this case the exponent in equation A.7T has a Valium 
ZeDG? COLE CS MONA Cao ai— mnie Actually @ in equation 
is replaced hy d, the spine diameter. Using thwvs value ge 

the general differential eyuation for the temperate 


excess bLecoues 
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Figure A.2 Cylindrical Spine 


{ 
where m = (4h/kdyt . 


ne general solution to equation A.4 is 


where the arbitrary constants are determined 


Heumdary conditions 


O(x=b) 
2 (= 0) 


eth 
‘DD 


(A.4) 


(A.5) 


from the 


(A.6a) 


(A .6b) 


This then leads directly to the equation for the tetperamume 


excess profile 


b cosh mb (A.7) 


and the heat flow entering the £1lnh (Spine) bases 


TT 2 


The value of h may be determined from equation (A.7) by 
taking temperature readings at the tip and the base which 
ace precisely "b" meters apart. 

Thus with base and tip temperature excesses, 

CiXSE =e elteeeay =a 


It is easy to see from equation (A.7) that 


a 


a cosh mb 
or that 


0 


aero Oley) = = 
a 
EEOUMmWIICh 
mb = arccosh(@, /@_) = (A.9) 
Then with a = (uh kaye one may then solve for h. 
Nae Gee 
ig ae 
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Figurce 2a) 


lcoy,Nu vs. logy,ka (0 degrees roll). 


32 


PO mee eee eee Be ete ee ee ee a ee a oe ee ee eee ee 


_ —- =» 
—- ae ee eee ew ee een ee ee ee ee ee ee en en ee —_— oe eee ee 
—_— 


LOGNU 
1.10 1.15 1.20 1.25 


1.05 


1.00 


0.5986R A**0.233 
ERIMENTAL DATA 


0.96 





meee eee eee eee ee ee — <r ee ee eee — <_ eee eees s ee e e ee e e e eee e 


2 lL eel 2 logaNu vse loy,Ra (10 degrees roll). 


os: 


—ee ee eee ee ee we wee eee ee ee ee eee ee eee ee ee eee eee ee ee ee ee ee ee eee ee eee oe 


LEGEND 
U = 0.6075RA**0.233 


N 
EXPERIMENTAL DATA 
RA 


e 


0/0 


LOG 





c2'T 0271 Ory Ort 90°F 00°! 96°0 


ee ee ee EE i ee ee ee ee ee ee ee ee ee ee ee oe ee ee eee ee ee —— << Fa a oP ee ee eT EE ee ee ee ee ee ey ee eee eee eee ee 





meg ee i eee ee ee ee ee ee ee es ee es eee ee es ee ee 


Figure 


leg, Nu VS. loj,Ra (20 degrees roll). 
34 


Urea 


—_ a SS a eee ee eee eee ee ee oe ee ee oe 


aa Se —— Ee eee ae eee eee ee ee oe 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
i 
{ 
| 
i 
| 
| 
t 


LOG,NU VS. LOG.RA 





Geol 


rae | 


O'T 
ON?9D07T 


907 


00'T 


EXPERIMENTAL DATA 


sp) 
+, 
nN 
© 
@ 
¢é 
< 
= 
Elen 
wy 
oe 
AS |S, 
I] 
= 
Fa 


O 
O 


$6°0 


LOG. RA 


— ee ee eee — oe — 
— -- ———_—_——— a ee wee ee ee ee ee ee ee eee ee eee ee eee ee ee ee 


Cure 


logaNu vs. log, Ra (Semucorees For). 


B.4 


S° 


Se ee Se 5 ee eee Se ee 
A SS ES ET EE Ay SE YY YE YE AS LY AS AY AY LY SS aS ET ee epee tae | 


JU2G RAO .2ae 


—_— eee 


MENTAL DATA 


fa 
ra 
x4 
OO 
Cx) 
rad 


LOG,NU VS. LOG.RA 
LOG.RA 





Get O2°T Qty Ory CO's 0o'T 96°0 


> ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee eee ee es es es ee ee i 


A ee ee ee <— ee ee wee ee eee ee ee ee ee eee ce ee ee ee ee ae eee ee 


(45 degrees roll). 


logpNu vs. lege Pa 


Bees 


Figure 


LEGEND 
0 DEGREES 
20 DEGREES 


C($) + LOG,RA 


Bm 
ty 
fe) 


0 


LOG,NU = 


- -_— —— wee ee ee me meee ee eee eee ee ee ee —_ om ear eee see oe ee eee ee i i ee ee 





oa eT ee 


owe ee ee —_ ———_ ST —Y 
a ee ee ee eee —e ow eee ee ee eee eee ee, eee ee ae ieee eee eet“ ee sm ee eee ee eee bad —— <> — — a —— oe eed 


Senprmrcon cf Deen Oemhesultse 
oI 


b 


Pace ee. 


7 
| 
| 
| 
| 
; 


(SaauDgaq) (9) ATIONV TION 





FIONV TIOU “SA ‘ON LTGSSON 


/ 
| 
i 
y OF OF ct oc C2 02 | fo | ¢ 0 | 
! . | 
' c 
| é : 
. Ces > WHY > 862g =a | 
QN399q91 ] 
\ | 
‘ 
| ! 
l ro | 
| Z ! 
' i 
Z. 7 
| y : 
YY) ' 
Y) i] 
se 
, no ! 
5 Pg ' 
) o O | 
| — ' 
> | 
J Cc | 
! = ! 
‘ ' 
H ‘ 
! I 
f 
: > | 
o 
i ] 
! | 
] 
l l 
; 
! | 
: | 
{ 
> po , 
3 
| 
{ 
] 
1 
| 


ve. Roll Angle (D) . 


Nu 


Eicjuc ee 


38 


geet x Ce 


TYPICAL DATA RUNS 


eid 


eines Cy. 
iypical Results (579 < Ra < 586) 
=—e790:.89'.% = 0° 
h 

4 oO oO 

n Tf K) ape @) Ra Nu NU oh (w fm? * °K) 
See Gh 2 8 ov det 0 2 eo 2 eo Gina 
3 ae Rp hare Ota) OF eel 2 sree Ze (eas: 
SA Zee U) See 979.40 Paes as, fae Si Chs 7 Oras 
Sy 2 cml Sis) oa amen oo Ose Zoe eu fore) 2 SAl 
S13 | Cyl 2 189s Seno Gg 2 wal Hele DAs) oa Shs WPS) (5) 8 
Cll”. Peoys) SSS Gia as 2.4842 24 So Ste) ois 
Cilia 208) Al Roles nec 2.4806 223005 6.843 
312.64 SieZae 4 20S aor 2./984 Zs om FE 8 
Ceara Siheaa luis 5 Orage 2.4066 2 fos el See) oir, 
oz CRI ails Ol) apes 2.7274 2. oon Tete 
SCZ ote Sole ek 3y 18 dls, elf 2 0 Sos 2.3304 142333 
Silay ders a Pats) Sys eae: Pago oat Vad is i ees 
See ee) Seleees 2 Seon as ¢ O94 Vai oye ess 7.446 
Beleve, Of Sills | oor Dara eas Olas Mage S) SIONS Gre 06 
BaleZ ic! Ses o 2 Oe Sl. Vier ch as 2.3340 fs 26 
Sy acme he) SB Ue abe 978.48 Pa es) Dono Peto 

meio 93 06CU C= CIP 
Sua 6, Seely, oS Oa zae 2 eee VANS io sine 7.934 
Blo Seo S) oalueele 2.9429 Ce Oro ed ieeOales 
Dee eer ak coal Bey 27) 2 OS 223g 96 oes) We 
Slee Slee anor Oe Jew 2233 Po SOS 
ol Zo Guar ole BiG 2 onl: 22072 PEBAC eae: 7.44) 
ool Zee Sore Src ish Gi 2.434 2.3406 ome) 7/0 
SZ ee Scll aors Oo Orcs 2.4434 2.3426 Gy. 7a 
3 2 26 Cn aces: 7oo suo ViMe lence: De aes Gs): 
Sa Belge ie 30270 Ze SOS Megs Sion! 67026 
Naa) cay S82 38 Se) gs oS IRS) F Cue 64 7.430 
SZ eG Se 7 Gots pol Saye emieil0) 3) 7 ieee a ay Yes al 
Sil Wat | ea ene Sie cure Zao. 23 tb 2.0 (BOs 
Sol? aie 2 Seo ao eyo Fee PR ASIS eps: Zoe HAG Tet 310 
ol ae SN ace oa pie 7 7 Ts eet ire ot Past 4 lois 
Slee, 6b Sa 2 7 ieee 2 2.17304 Tapas ote, Ts Sree 
See Y1 Ce S84 .44 eer. 3 0) Des orth 74004 


dyes gle 


Cra 


Typical Results (579 < Ra < 


= 296. 


Soler. 
oi eae 
Sols 
Slee 


3 eee 


joe 
ols 
ieee 
eile 
3 ae 
Suze 
Se 
Cae 
Sel Dis 
cWllyae 
Sua 


084, 
5) ee 


sien: 
Siefisi 


See 
S51 oh 
Show. F 
Dect! 
Solem 
of Ele. 


2 Oe 


Bile tour 
eee) c 
OBee 
Jor 
Shot leg 


Some 


40 


POND ROCIO ORD RO oho ts) PD ei iat 


Soho Uhoobor BS ID ROO ND RO RD IND NOOR ND ho 


586)) (eenmuaaaee 


7 Sed 
. Oho 
» 2 2aNb 
» food 
Ree Oke 7. 
Bele los 
Hi celell 
Be sus 
Be chiel 2) IL 
:odleew) 
~6246 
Bas) ee) 
p63 
. 7680 
~7584 
Pete sy ci 


~ 7480 
oe 
Boye lr 
. 7684 
slione 
Bbc he ede 
Yor 
oe 
Pe lkss6 
2 oa 
~6 842 
Bots) s) 
Ory 353 
rie cous i 
Me ee! Ah 
Mipsvcnl 


NOUNS ROU Roe fh BRO hOr hho th ho ahaa 


Sere hh Oh > I BO RO ND RD NM) NO RO BO 


sou 
3332 
. o Som 
«3325 
-ocus 
- 3304 
soe 
3 eo 
os 4s 
vo ous 
Merch s) |) 
. 3406 
Pees S| | 
 o828 
~3 oe 
23345 


. soa 
: 3c 
Moos 
. 3 oa 
. 333i 
~3384 
3 3oe 
. done 
3 oe 
3 3am 
23303 
3408 
32a 
32a 
~3 344 
~3 344 


Sn nNnynnr orn ononrnnO~) ~ 
se oe @® #@  @ ¢- ee @  @®  @ ee e@#  oe® @ 


SYA NYNNANMDANIDANW NAIM) ~I 


- Sa 


. 608 





wos Oo ££ we EH 


Oe: 


See 1 


iapuead ire cmure (579 < Ra < 586) Ceont'd.) 


Yl 


IND IND ROOD NOOR DOTS NO NSD re oh 


~21640 
oe, Paar 
Bo) cool a 
2 
7a SM 
~4b9s 5 
~46626 
22S OS 
so S.C 
et ae OS) 
sibig oe) 
Hs oy SIO, 
(VSO 
2h OG 
ies ee Oy 
. 76402 


INS ND NSN RD Roh NS Os IS PS 


cis) ors) 0) 
fee) Com 
7 oS OIG 
Meio 0) FAL; 
Sen ios 
Mpcks 6 Silo 
Reo orolisns 
 SpopenleZ 
ore 0 
0 07 Oo 
eo lea 
Peas 
poe / 
ieee D 
A005 6 
Ho Oro 


JIn~nInNnNInNANOMOOOD IJ I~) ~) MD 


eS 
ye ou 
. 004 
.o081 
owl ual 
Bie) lee) 
nou 
meee 4 
~594 
ore |. 
moo 
“076 
eS 
./709 
Mees 
v0 25 


gpa gle |e doe Cca et 


Typical Results (1305 < "Ra <cee 


Jt ee ee aces, * ss 
S 
S @) 
Pin Ty, ¢ K) 
Hi 3369 
2 338.49 
3 3S Gro 
y 338.44 
5 336 oss 
b gee 
7 Scr 
8 3302250 
g 366.00 
130) Spots emer. 
ama Bie cm ao my 
V7 33 eeGr 
13 3 3'8-25.0 
Vy So oreo: 
5 Sia cnes ts 
ike 33640) 
dee =) ue 
S 
il SoHo 2 
2 See Ba 
3 See eet onal 
uy eretae os. 
J 34a 6 
(, Soden 2 
7 S Siiteg © 
& Be Pasiis! 
0 3 Stee 
10 Sore. 
a: Soa 
eZ SisAbes Cee 
ine cel SES 
14 Gee eertel @ 
ee 334253 
ive Cee 30 


Bones |e! 
By eus) 10s 
33 Geez 
Cc Al 7 
33320 
Rpome (0)! 
330.00 
Spence Olt 


LG? 


334835 
Be hey, 
S34203.3 
3 3s 5 
Sos 
33H ss 
33 0225 
Ge less 
33 heh 
3342.37 
334.4] 
Ghee ol 
334 2.49 
Seen 
33 teen 
33 lias 


Ra 


1304. 
30s 
sition, 
JE SUenO 
nS OLS, 
pesi(0) 5 
PsUse 
dies 0c + 
SiG Ga 
sree. 
LSaeee. 
alms elie 
lotsa 
1304, 
1 S0ice 
ES Uae. 


Lo Zi 
Lsaee 
Ne 2 
Loge 
L323 
Loe 
Looe 
fe 5) 
hase 
lS Zee 
Ls 23e 
Ue) cil, 
Teo2ue 
1 32 on 
3 Gee 
17 oe 


42 


GO G) GO G) G@) NO NM GG GW @® GW Wt Gro 


WW WWW WW WWW WWW WW W 


Nu 


374s 
20 oe 
.07919 
~28048 
.15390 
. Lares 
.11699 
. 2O2 one 
meee ois 2 
.96674 
98347 
.06200 
\ 2 oa 
-O45 04 
 beeeD 
SoU Gt 


2 oe 2 
Umer 2 
2 Ole 
<0 43a 
i OE 
00764 
oe Ge 
3 Sed 
sO 3803 
s0 56:39 
» 2.0 Oe 
eres 17) 
: ibss i 
«see eee 
eZ 
S20 508 


NO RONDO RO ROE RON NOON TO NOES oh hor 


NOI NI RO RDN MH DN NM NN TIN MN NN 


Nuon 
. (OU Sis 
.76029 
. os 
. (550s 
76146 
(oO se 
~. oles 
.76064 
. 7 Omer 
, -omree 
. Posey 
76420 
.76064 
. oie 
. /Gleas 
, FOU 


«7 Toe 
1 1254 
«7 1 Se 
Rey blo) / 
. (Oo Sree 
i 7 ee 
. (6 Sia 
1 i 2m 
7 12 
of Flee) 
712s 
- 1/1302 
7 Peo 
/ Oa 
Mees! 
71 Oge 


(w/m- ee 


>on 
2 
~78Y 
22m 
~949 
. o oe 
. oom 
. L224 
. 6 ike 
. 5am 
Moy 
. Pee 
. Lz 
. i 
mee 5 
2 oe 


JITIONDOODDOD YH MH OD ~~) O ~) O 


WAMOMDANMAMMWAIAIAANDODODNANANM ~) 


h 


. 008 


Typical Results (1305 < 


ie S11, 24 
S 


B50, 
SS1Ie 
She 


 e= 


IO 


Sou Oe 
So 0e 
3303 
Cee ce 
330z 
ems (Ul 
3302 
i611 9)7) 
ge. 8) 
o7r 
Soe 
Sic Ue 
53.0%. 
550. 
S30e 
S72 oe 


TAS be Cc. 2 


au ee 
1340. 
bie Sie 
P3sb% 
es eie 
ee 
HDS os 7 ee 
eae 
1s 2 Sie 
Ne sila 
is ow 
1346: 
1342. 
Sco 
see 
ror 


eS ae 
eS Ie: 
Lc ae 
eS aor 
ee) 
eo Ze, 
Nei 2 
IRs) 7s) 
fies onli 
NG Siem 
ye Gs: 
blige 22'S 
iwc 
cy oe 
Nee soe 
i vec 


43 


hee 2 Geen tk.) 


OCI '°GO- Go CO Go Go Ol CoCo Go Go Oo Ged Or 


fa Oo Co CO ECO CoO COG Oo CO) Co Co IN 


~11643 
7 QI 
- 2 tie Oro 
raver lite 
2 oe 
eo OS a7 
Rae) Slice 
ele oie hs) 
~20ib 3 
V2 Zot 5 
23S 
Ape wis) chs 
rare 0) Jeus) 
fay he 3) 
Geese 
Ba ca!) 


elope hin 
Reel 3) chile 
mlbic Oise 
lc Oring 
plies Silene: 
me elon 
~YYUIGY 
Bd oy ayaa 
vo oes 4 
vie oe 
lke wee) 
ele) chs V7 
4S 36:0 
oo 29 
Sew 0 
5 dea Mes 


hI NO ND ND RO TROD RO NOON RD SS Re hh 


Moh OND RO RD TORO ND PO RO Nh RO PS 9 


Pe se Oe, 
77741 
.77704 
Me ioucys. 
ae 23) 5) 
Me itera! 
A esas 
per es 
Be iiscons 
Tis Se 
nat OT 
eS 
sa) Wecne2 
oe Wand 
<7 SO 2 
i eels. 


oi S139 
Bes) ots 9 
a oles 
ce SeleeS 
reece line 
me aoa. 
eorise I 
5 STEN OPS 
pea del 
oa 
19514 
mle elleyr4 
Rye cleo ts 
ps2 7 
Seco ke 
eo aaa 


MOAYDIDMNOMTITAIAINAINDMMDNMD MOM 


TOMO DAINIOIOOIOOANANDMD MM 


eG 


628 
453 
823 
823 
887 
082 
yea 
uae 
625 
.799 
760 
.845 
S44 
669 
ite 
aye 


TARBLE Gee 


Typical Results % 03 qo 


c= oO SO ne 
S 


Soe, 
S 


3 eo ecoes 
3 oe 
oe 

cus ies le 


deve bee ae 
eo rer 
i Sales 
digs sei 
Usa Ee 
IL ade 
eee oe 
dbs) Jeep 
esol: 
dete ne 
AS rales a 
lees 
eS 2s 
eS 8e: 
pipes evans 
IRS. 


Le3e2 Uae 
gs 204 
ie ee eee 
gies 
1322 
eo ee 
seca es). 
Vo2ee 
ils 
L313 
32 ee 
LS e 
is ibe 
See 
SS 
L320. 


yy 


GO) GC) GN GN ® DW GO GON (GD Co co ce 


ee CA Gade nea) GOI IND CO) (G0 GO) Gad (Ca Ca) CG 


< i380) 


»3Pase 
Oe Ghs 
»29894 
id vee oe 
~16447 
Me ietos:) fo 
a oe 
- 20 C02 
Bc 70 
eee Sal 
sO 6 
sO eter era 
pales Ih? 
se LES 
42424 
stoma 


~10591 
, 10 ae 
SAG 
C4478 
JONES 2 
<0 64 dei 
Piles 
Okc Ron Ons) 
Pee o) Shes 
~9540s 
,Ou56u 
429 18 
~48656 
3 esas 
prs sng 
CORES 


ND ND ROO ND BOORD NORD ROR NO DOs oo 


oe ere IND TS OD Dh RRND RD ROD ND 


Ceconte. 


. 1S ee 
(oes 
: /oeeial 
. | OSD 
~-76824 
. / Cee 
: 1 6iGaw 
. ote 
. 76806 
. (Obra 
. Cae 
- } Gai 
. (oS 34 
.) ees 
, Poa a 
. / Gdaes 


. (6Bee 
. 69m 
7 1 Ome 
1 TOs 
76989 
,  orOs 
, ogi 
Rasy?) S| 3) 
. | Pears 
. (oa 
<1 Feo 
, (Br eZ 
. boOsG 
, 6 ial 
, | boa 
. 76904 


~~) O&O Ooo &) & OF Oo 12032 2) 


ead) <) ~) OF OC? OC) > &) © O&O) O) O& O Oo 


See 
~468 
.864 
. Gis 
. Gigns 
2 oe 
. Dae 
. Siew 
. See 
. Saez 
. 2o 
~400 
~246 
. oe 
. 2a 
44g 


462 
~4YO4 
~664 
. Sees 
~344 
. Saas 
~466 
448 
. om 
~14E 
mec) 
. Loe 
2oe 
. 02s 
_ oe 
~4BY 





Aer ©. 2 


ieprcadebeomeus (lots < Ra < 1380) (Cont.) 
T.=301.86 i = 20,05 
7} gue oul 318.58 ieee PSB ee 2 7 9NES 5.869 
2 318.87 ese ieee 3227655 62. 79204 6.884 
3 1396 318.56 Mev ome 8) ese oo 65o ~ 2. 7/9480 pen) OF 
r 318.96 318.57 1379.99 3.21071 2.79489 6.747 
C, 2 eh le 318.62 Megs)? Soev0213 2.79627 a BS 
6 318.97 318.59 eciero 8 3.70759 2.795450 6.740 
7 2 Ee Oy, S168 lcci tomo e4so5). 2.79669 6.823 
g 318.95 318.56 ieee oc geeroe 7 scss | 6279462 6.850 
g 318.89 318.52 1375.58 3.12968 2.79298 6.576 
10 378.94 318.57 ic ee eo lO) )2a794 71 6.554 
1] 318.95 Sees eyo Sewer eo ) 2h 7 1 6.749 
1 2) Tate eo Ma 318.58 eco so Oona 2S 83 7.426 
13 318.99 oe. 55 Were 6. 8.58582 2.79572 5 BS 
14 318.96 318.59 ecm domes, 02125. 2.7952? 6.349 
15 Slee: cu 318.54 ies) Mees 07 7 ome 2 ae 37 7 6.466 
16 318.98 318.59 NCS. Some ecm 2 oS 4 8 6.838 
T,=309.50 i ee 
] oO ae mall S779 eee 7 o Wes O7 oll) 2.77703 6.6)] 
g 328.26 B27. 83 Penge eh wee O87 2 9.77833 6.866 
3 a7 Ree | 32ers 7 ess 2 oi 9 7 887 6.406 
ly 328.16 ee 2 ie cocen so asco le? 2477532 eS 
4 eS B75 Meccr ot Ce OdGn 8s @257 7562 6.450 
7 CUES ere B27 484 Motmueeg pee gn 55 9.477772 5.793 
7 20 ool 327.84 otmeesGe eee 899 9.77772 5.974 
g 878231 327.93 oem o OUTS. 2-7 8042 6.029 
G one tl S27 = Fal Mie eee? 8 0 77 58S 7.980 
10 a28 . 09 227 Jak iPoeeeoo) Je eol0S §2.77431 6.193 
at eee 6 ee So ousioess cetbiysl 9 2.7777? oS ae: 
1? 328.76 Goes t he cles sce 46963 2277780 7.502 
ils 37 Bye 227.63 eee). se 75752. 2.977547 e207 
1y 32 2a 627.72 Neca oo e132. 2.77469 6.279 
‘5 328.26 Oy ere eee? 59 93527345565 27277826 6.95 
16 be 8. 21 S27 S ieoces Seo euSGgs 22977648 Fone 


4S 


(een ey 0) CC) - CCC) s-Co C2 C2 Co OO OO Oe a OC a OC Cr. C2 C0 02a aa 


APPENDIX D 


DATA REDUCTION PROGRAM 


ENE NISES By ee aera eee re Aree ars OF ai6 OK AA re AYE AS es Bs OS OES 8 M4 MMM ™M MM 
WWHRWRWRW WO ORO KRRRKERRREFRREKKKHKKRHEMHHKEKHI THHHKHKKHKHHKHHHKH 


PROGRAM TO CALCULATE HEAT TRANSFER COEFFICIENT ON A : 
ieee Eee DISK % 


a 


¢ Ae * Ae * Ae * ee | ee | 2 Y 2, YAY A, © A, Y 2 YY A YY 4 fs 
MHKHEKHGEHMHHHKHHHKKHHKMX MORK HH STH HRS HHH HRHMHHEMHHRHRHKHKHKHRKHKHKRKKXK 
f aes, © vi 
RREKKKKEKRKRRRKKKKEDRET NITION OF VARTA BLE GARR KRKRKKRKKKEKRKKKKRES 
pA = FACTOR IN CHURCHILL AND CHU CORRELATION USEDSRas 


TO SIMPLIFY AND SHORTEN THE EXPRESSION 


B = LENGTH OF THE SPINE 

C = SIMPLIFYING FACTOR FOR CHURCHILL AND CHU CORRELA- 
TION 

CB = CORRECTION TO THE SPINE BASE TEMPERATURE 

GO = CORRECTION TO THE SPINE TIP TEMPERATURE 

D =- DIAMETER OF SPINE 

E = SIMPLIFYING FACTOR FOR THE CHURCHILL AND CHU COR- 
RELATION 

¢ = GRAVITATIONAL ACCELERATION 

GR = GRASHOF NUMBER 

H = HEAT TRANSFER COEFFICIENT (EXPERIMENTAL) 

KA = THERMAL CONDUCTIVITY OF THE FLUID (AIR) 

KP = THERMAL CONDUCTIVITY OF THE SPINE MATERIAL 

M = A HEAT TRANSFER PROPERTY 

NU = NUSSELT NUMBER (FROM EXPERIMENTAL DATA) 

Wc ~ = a NUMBER (FROM CHURCHILL AND CHU CORRELA- 
TION 

NUAVG = AVERAGE NUSSELT NUMBER (FROM EXPERIMENTAL DATA) 

NUT = TOTAL OF NUSSELT NUMBERS FROM 16 SPINES (1RUN) 

PHI = ROLL ANGLE (DEGREES) 

PR = PRANDTL NUMBER (EXPERIMENTAL) 

R = RATIO OF TEMPERATURE EXCESSES ( = THB/THO) 

RA = RAYLEIGH NUMBER (EXPERIMENTAL) 

RAT = TOTAL RAYLEIGH NUMBER FOR 16 SPINES (1 RUN) 

RAV, = AVERAGE RAYLEIGH NUMBER (EXPERIMENTAL) 

TB = BASE TEMPERATURE (DEGREES F) UNCORRECTED 

TO = TIP TEMPERATURE (DEGREES F) UNCORRECTED 

TBF = CORRECTED BASE TEMP. (DEGREES F) 

TOF = CORRECTED TIP TEMP. (DEGREES F) 

TB(I) = CORRECTED BASE TEMPERATURE (KELVIN) 

TO(I) = CORRECTED TIP TEMPERATURE (KELVIN) 

THB = TEMPERATURE EXCESS AT THE BASE OF SPINE 

THO = TEMPERATURE EXCESS AT THE TIP OF SPINE 
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TEMPERATURE OF FILM BETWEEN SPINE AND FLUID 
TEMPERATURE OF SURROUNDING FLUID (AIR) 


*DECLARATION OF VARIABLES* TRS 


REAL D,KP,B,TS, Pee eee NU Te nO. Rael, RAV, ‘NUT, ‘NUAVG, 


Ri A Cenc ln VA, ALFA, KA, GR,PR,RA,NUC 
REAL TB(2 Dy hee), 25a G0(25) 
INTEGER I 

SE 3 3 SE Cree Ni | Peek | LON OF VARIABLE SRR RKKEKRREEKRER 
BS 020127 
B = 0.0508 
KP = 120 
ee 9.807 
To S753 
PHI = 45 
Mo =20 
RAT = O 

CALCULATE TS IN KELVIN 
moos gS = 32)°5/9 + 273.15 

PRINT HEADINGS 
WRITE(6,100) TS, PHI 
WRITE(6,150) 

READ IN AND CALCULATE CORRECTED TEMPERATURES, DEGREES F 
DO 50 I = 1,16 

Raise eOe) TB Cr), TO(T),CB(L),cOo(1) 
TBF = TB(I) + CB(T) 
cOre= TO) + COLT) 


DONVER CORRECTED BASE @ TIP TEMP. TO KELVIN 


TB(L) = (TBF - 32)#5/9 4273.15 
moO Tom TOR = 32)%5/9 4273515 
CALCULATE THETA BASE AND THETA ZERO (TEMP. EXCESSES) 
Rio eT) = Ss 
Migt=.TOlT) = TS 
Pvc UMATE RATIO OF THETA B TO THETA O (=cosh(MB)) 
R = THB/THO 
CALCULATE M 


M = (ALOG(R + SQRT((R**2)-1)))/B 
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ee ee ap ee 


CALCULATE HEAT TRANSFER COEFFICIENT H 
H = (M#*2)*KP*D/4 
CALCULATE FILM TEMPERATURE 
TF = (1./4.)*(TB(L) + TO(L)) Gin 2a eae 


CALCULATE PROPERTIZS SORA ais 


KA = ©.004033 45072 20t- us 
VA = (-1.22567 + (9. 350-3) *7R) = caeee 
ALFA = (-1.82333 + (0.01350°TH) (7 Gem 


CALCULATE NUSSELT NUMBER BASED ON EXPERIMENTAL H 


NUS = aby 
NUT = NUT + NU 


CALCULATE GRASHOF, PRANDTL, RAYLEIGH, AND NUSSELT NUM- 
BERS BASED ON CHURCHILL & CHU CORRELATION 


PR = VA/ALFA 

GR = (G*(1.0/TF)*(2P-TS )*D**2 7 

RA = CR 

A = 0.387*RA**(1.0/6.0) 

C= (1.0 + (0.559/PR)**(9.0/16.0))**(8.0/27.0) 
fe = he 


NUC = (0260 2) 2 
Rad eRe tee 
PRL Oa wal 


WRITE(6,250)1,TB(L),1TO( 1) yee Gee 
CONTINUE 


RAV = RAT/16.0 
NUAVG = NUT/16.0 
WRITE(6,300)RAV,NUAVG 


STOP 
FORMAT('1!,15X, 'TS=',F6.2, 3X, 'PH1=',F4.1) 

FORMAT(17X, 'PIN', 4X, 'TB', 7X, '70', 7x y RA, Coogee 

4 'NUC', 7X,'H') 

FORMAT(4(F 23)) 

FORMAT (18X,12,2X, F6. 2, 3X, FG. 2, 2X, F7.2, 2; Poe5 ye eee 
# 3X, F6.3) 

FORMAT(16X, 'RAV =',F7.2,2X,'NUAVG =',F6.3) 


BND 
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